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摘 se. 喀 喇 昆仑 公路 是 推动 我 国 新 疆 地 区 对 外 开放 的 重要 通道 ,长 期 受到 冰川 灾害 的 影响 , 尤 以 罕 萨 河中 游 的 
巴 托 拉 段 最 为 严重 。 本 研究 基于 野外 考察 .和 感 影像 等 资料 ,分析 了 巴 托 拉 、 帕 苏 ` 固 尔 金 冰 川 近 百 年 的 变化 .变化 
原因 及 对 公路 的 影响 。 结 果 表明 :(1) 近 百 年 来 ,特别 是 2000 年 以 来 , 巴 托 拉 、 帕 苏 、 固 尔 金 冰川 总 体 处 于 微弱 退缩 


状态 ,对 公路 的 直接 威胁 趋 于 减弱 。(2) 巴 托 拉 、 帕 苏 、 固 尔 


金 冰 川 消融 加 剧 与 冰 舌 波动 导致 融 水 径流 频繁 迁移 , 冰 


川 融 水 洪水 ` 冰 湖 溃 决 洪水 .冰川 泥石流 等 次 生 灾害 对 公路 安全 运行 的 影响 正 日 渐 凸 显 。(3) 巴 托 拉 段 冰川 变化 对 
公路 最 直接 的 威胁 在 于 2021 年 6 一 7 月 巴 托 拉 河 新 的 改道 ,改道 后 水 流 将 直接 冲击 侵蚀 路 基 ,影响 公路 安全 。 本 研 
究 成 果 可 为 中 巴 公路 治 线 以 及 我 国 西部 与 中 亚 .南亚 其 他 陆路 通道 的 冰川 灾害 监测 防治 提供 参考 借鉴 。 
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喀 喇 昆仑 公路 (Karakoram Highway ,KKH ,以 下 
简称 公路 ) 北 起 新 疆 喀 什 , 南 至 巴基斯坦 塔 科 特 , 是 
中 巴 两 国 实现 互联 互通 ,建设 中 巴 经 济 走 廊 的 战略 
通道 。 公 路 穿越 的 东 帕 米尔 高 原 . 喀 喇 昆仑 山 是 高 
亚洲 冰川 分 布 最 集中 的 地 区 之 一 ,冰川 灾害 曾 多 次 
威胁 和 中 断 公路 通行 ,造成 重大 损失 "”。 位 于 巴 基 
斯 坦 北部 罕 萨 河 ( 洪 扎 河 ) 中 游 的 巴 托 拉 段 是 整个 
公路 沿线 距离 冰川 最 近 、 冰 川 灾害 最 集中 的 路 段 ， 
其 中 以 巴 托 拉 (Batura) 、 帕 苏 (Pasu) 、 固 尔 金 
(Ghulkin) 等 3 条 冰川 对 公路 的 威胁 最 大 也 最 为 直 
pe^, 

1973 年 , 巴 托 拉 冰 川 主 排水 道 ( 巴 托 拉 河 ) 改 道 
引发 的 冰川 洪水 冲 毁 了 建设 中 的 公路 和 公路 桥 。 
1974—1975 年 ,为 提供 合理 的 修复 设计 ,中 国 考察 
组 对 巴 托 拉 冰 川 开 展 了 包括 末端 进退 流速 .厚度 、 
融 水 径流 等 全 方面 .详尽 的 科学 考察 ”。 最 终 考 察 
组 提出 按照 原 线 修复 的 科学 建议 ,确保 了 公路 的 长 
期 稳定 通行 ”。 此 后 的 1978 年 1980 4E. , 1994 年 
等 ,考察 组 再 次 对 巴 托 拉 地 区 开展 科学 考察 ,积累 
了 丰富 的 地 质地 貌 冰川 物理 .水文 .气象 等 科学 资 
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料 "。 在 此 基础 上 ,学 者 们 揭示 了 巴 托 拉 帕 苏 L8] 
尔 金 等 冰川 在 20 世 纪 近 百年 的 进退 变化 ,并 预测 巴 
托 拉 冰 川 将 于 2010s 前 后 再 次 发 生前 进 “ 7 。 

进入 21 世 纪 以 来 , 随 着 气候 波动 的 加 剧 , 巴 托 
拉 、 帕 苏 、 固 尔 金 等 冰川 发 生 了 许多 新 的 变化 ,这 些 
变化 对 公路 的 正常 通行 产生 了 新 的 威胁 。 例 如 巴 
托 拉 冰 川 的 主 排水 洞 位 置 不 断 北 移 ,最终 导致 融 水 
通道 被 冰 舌 撞 奈 切断 ,时 隔 近 50 a 后 再 次 引发 巴 托 
拉 河 的 大 改道 ; 帕 苏 冰川 末端 的 冰 夺 阻塞 湖 频繁 演 
决 ,2008 年 1 月 和 4 月 2 次 引发 冰 湖 演 决 洪水 ,中 断 
公路 通行 ; 固 尔 金 冰 川 末 端的 频繁 波动 导致 冰 崖 
下 的 排水 通道 不 断 迁移 , 仅 在 2008 年 就 引发 4 次 泥 
石 流 ,致使 公路 桥 底 严重 汶 积 “小 。 因 此 ,鉴于 灾害 
风险 的 日 渐 凸 显 和 公路 无 可 取代 的 地 位 ,有 必要 对 
巴 托 拉 、 帕 苏 、 固 尔 金 冰川 的 历史 与 新 近 变化 进行 揭 
示 和 分 析 ,为 公路 的 安全 运行 和 灾害 防治 提供 参考 。 


1 研究 区 概况 
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西 段 主 疹 一 一 洪 扎 喀 喇 昆仑 山北 坡 ( 图 1)。 洪 扎 喀 
喇 昆仑 山 平均 海拔 超过 6000 m, 最 高 峰 巴 托 拉 峰 海 
拔高 达 7795 m, 为 冰川 发 育 提 供 了 极为 有 利 的 地 形 
条 件 ,使 得 巴 托 拉 冰 川 群 成 为 了 公路 沿线 规模 最 
大 、 末 端 海 拔 最 低 、 进 退 最 为 频繁 的 冰川 群 之 一 。 
巴 托 拉 、 帕 苏 、 固 尔 金 冰川 上 游 粒 雪 伟 面积 开阔 , 补 
给 物质 充足 ,中 下 游 则 深入 了 干 热 河 谷 , 消 融 强 烈 ， 
因此 属于 典型 的 多 温 型 冰川 ,并 表现 出 较 高 的 物质 
平衡 水 平 “。 本 研究 的 解 译 结果 表明 ,2021 年 巴 托 
拉 冰 川 长 约 56.5 km, 面积 307.12+4.32 km? ,未 端 海 
IR 25 2620 m; 帕 苏 冰 川 长 约 26.4 km, 面积 64.81+ 
0.72 kn? ,未 端 海拔 约 2700 m; 固 尔 金 冰川 长 约 17.8 
km ,面积 28.26+0.47 km ,末端 海拔 约 2540 m. 


2 数据 与 方法 


2.1 数据 源 

2.1.1 Landsat 影像 ”本 研究 使 用 了 获取 于 1972 一 
2022 年 的 25 景 Landsat 系 列 影像 ,以 揭示 研究 冰川 
的 相关 变化 ( 表 1)。 所 用 影像 均 为 LIT 级 产品 ,由 
USGS (United States Geological Survey ,https:Wearthex- 
plorer.usgs.gov/) 进行 了 系统 辐射 校正 和 几何 校正 ， 
并 结合 DEM 进行 了 地 形 校正 。 由 于 Landsat 系 列 影 
像 有 相当 高 的 正 射 校正 精度 ,达到 1/2 像 元 左右 ， 
因此 本 研究 直接 利用 Landsat 影 像 开 展 冰川 变化 研 
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究 ”。 为 提高 冰川 边界 .出 水 洞 . 融 水 径流 等 地 物 
的 解 译 或 识别 精度 ,本 研究 选取 了 消融 季 末 期 ok 
川 区 无 云 .山体 阴影 覆盖 较 小 的 影像 ”。 

2.1.2 冰 面 高 程 变化 数据 集 Brun SE AJH ASTER 
立体 像 对 提取 DEM, 基于 大 地 测量 学 方法 ,通过 
DEM 的 空间 匹配 与 误差 校正 ,重建 了 2000 一 2016 年 
高 亚洲 冰川 的 表面 高 程 变化 与 物质 平衡 ,并 生成 了 
高 亚洲 冰 面 高 程 变化 数据 集 。 数 据 集中 的 误差 主 
要 源 于 DEM 的 生成 精度 及 其 空间 匹配 偏差 ,因此 以 
非 冰川 区 等 地 形 稳定 区 的 高 程 残 差 来 校正 并 评价 
冰 面 高 程 变化 的 误差 “"。 该 数据 集 的 空间 分 辨 率 
为 30 m, 可 获取 于 PANGAEA (https://doi.pangaea.de/ 
10.1594/PANGAEA.876545)。 

2.1.3 冰 面 流速 数据 集 NASA 的 ITS_LIVE (Inter- 
mission time series of land ice velocity and elevation ) 
产品 包含 了 1985—2020 年 高 亚洲 主要 冰川 的 流速 
数据 。 该 产品 基于 Landsat 影像 ,采用 Gardner 等 提 
出 的 自动 裂缝 特征 跟踪 处 理 链 方法 ,通过 进行 局 部 
归 一 化 .过 采样 ,特征 跟踪 等 提取 冰川 流速 *”“"。 影 
像 的 空间 匹配 偏差 是 该 产品 的 主要 误差 来 源 ,尽管 
在 生成 过 程 中 已 基于 基 岩 等 稳定 地 形 对 其 进行 校 
TE ,但 无 法 彻底 消除 。 因 此 ,ITS_LIVE 产 品 同样 基于 
非 冰 川 区 的 流速 残 差 评价 冰川 流速 的 不 确定 性 。 
1999 年 之 前 (CTM 影 像 ) 冰 川 流 速 估 算 的 不 确定 性 超 
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图 1 巴 托 拉 、 帕 苏 、 固 尔 金 冰川 的 地 理 位 置 
Fig.1 Geographic location of the Batura, Pasu, and Ghulkin Glacier 
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表 1 遥感 影像 数据 
Tab.1 List of the Landsat images used in this study 
影像 编号 获取 日 期 /年 -月 -日 分辩 率 /m 影像 编号 获取 日 期 /年 -月 -日 ”分 状 率 /m 
LM21610351972251AAA05 1972-09-07 60 LC81490352013282LGNO0 2013-10-09 15 
LM21610351979250A AA05 1979-09-07 60 LC81500352014276LGNOO 2014-10-03 15 
LT51490351990219ISP00 1990-08-07 30 LC81500352015231LGN00 2015-08-19 15 
LT51490351993163ISP00 1993-06-12 30 LC81490352016275KHC00 2016-10-01 15 
LT51490351994246ISP00 1994-09-03 30 LC81490352017261LGNOO 2017-09-18 15 
LT51490351997222SGI00 1997-08-10 30 LC81490352018216LGN00 2018-08-04 15 
LT51490351998225XXX01 1998-08-13 30 LC81490352019267LGNOO 2019-09-24 15 
LE71490352000255SGS00 2000-09-11 15 LC81490352020238LGN00 2020-08-25 15 
LE714903520012418GS00 2001-08-29 15 LC71490352021160LGN00 2021-06-09 15 
LE71490352002276SGS00 2002-10-03 T5 LC71490352021192LGN00 2021-07-11 15 
LT51490352008285KHCO1 2008-10-11 30 LC81490352021272LGN00 2021-09-29 15 
LT51490352009223KHCO00 2009-08-11 30 LC81490352022147LGN00 2022-05-27 15 
LT51490352010290KHCO00 2010-10-17 30 
过 10 ma',1999 年 之 后 (ETM+ 影 像 ) 降 低 为 约 2m: 的 像 元 数量 ;A 为 像 元 面积 。 
a',2013 年 之 后 (0LI 影 像 ) 降 低 为 约 0.5ma 7*7, 2.2.3 出 水 洞 与 融 水 通道 识别 ”出 水 洞 指 位 于 冰川 


本 研究 选用 了 其 中 1985—2018 年 的 冰川 年 平均 流 
速 合成 数据 ,空间 分 辨 率 为 120 m (hitps://its-live.jpl. 
nasa.gov/) o 
2.2 方法 
22. 冰川 边界 解 译 ”由 于 研究 冰川 表面 存在 大 范 
围 表 研 覆 盖 , 且 进退 变化 频繁 ,本 研究 采用 目 视 解 
译 提取 冰川 边界 ”。 在 中 上 游 裸 冰 区 ,直接 根据 
Landsat 假 彩 色 影 像 ( 热 红 外 、 红 、 绿 波段 组 合 ) 中 显 
著 的 色彩 差异 进行 冰川 边界 的 数字 化 (图 1)™。 
在 表 磺 覆 盖 的 冰川 下 游 , 则 根据 末端 冰 岩 与 出 水 洞 
的 组 合 、 高 程 变化 等 地 形 和 水 文 特征 识别 冰川 边 
FR?) BOR ae RP ROK KG Fem KK 
EAK 、 冰 崖 与 水 流 等 的 色彩 差异 均 得 到 了 显著 增 
强 ,由 此 可 准确 识别 冰 舌 、 出 水 洞 等 的 位 置 。 野 外 考 
察 资料 表明 ,研究 冰川 冰 舌 消融 强烈 ,因此 末端 冰 岩 
与 出 水 口 的 组 合 即 代表 了 冰川 运动 的 最 前 缘 “”。 
在 实际 解 译 中 ,本 研究 首先 完成 2021 年 的 冰川 边界 
提取 。 在 此 基础 上 结合 不 同时 期 的 影像 ,修改 冰 舌 
等 发 生变 化 的 部 位 ,以 获取 其 他 年 份 的 冰川 边界 。 
2.2.2 解 译 精度 评价 ”基于 光学 影像 解 译 冰川 边界 
的 精度 主要 受 影像 空间 分 辩 率 的 影响 ,因此 本 研究 
采用 统计 冰川 轮廓 线 经 过 的 像 元 数量 的 方法 ,评价 
目 视 解 译 提取 冰川 边界 的 不 确定 性 ””。 公 式 如 下 : 
E=N-A/2 
SUH: E HIKI ES AVR RE TE 5 NK Gd Fe 


末端 , 冰 内 和 冰 下 水 系 的 排水 出 口 , 本 人 研究 结合 野 
外 考察 记录 和 遥感 影像 识别 冰川 出 水 洞 与 融 水 通 
道 的 迁移 变化 。 如 图 2 所 示 ,在 消融 季 末 期 的 Land- 
sat 和 Google Earth 影像 中 ,得 益 于 巴 托 拉 冰 川 较 大 
的 排水 量 , 冰 川 出 水 洞 与 融 水 通道 的 准确 位 置 可 获 
得 有 效 识别 。 此 外 ,在 1973 一 2021 年 2 次 大 改道 期 
间 , 巴 托 拉 河 下 游 保持 稳定 ,上 游 随 冰 舌 波动 而 延 
伸 , 因 此 结合 不 同时 期 的 遥感 影像 即 可 准确 识别 出 
水 洞 与 融 水 通道 的 变迁 过 程 。 


3 结果 与 分 析 


3.1 冰川 进退 变化 

近 150 a 来 , 巴 托 拉 冰 川 发 生 过 多 次 进退 波动 。 
1873 年 巴 托 拉 冰 川 发 生前 进 ,但 此 后 近 50 a [E] os 
位 置 基本 保持 稳定 。1930s 冰川 开始 衰退 ,1944 一 
1966 年 退缩 827 m。1960s 末 转 入 前 进 ,1966 一 1980 
年 前 进 约 136 m。1980s 初 再 次 衰退 ,至 1994 年 累计 
后 退 了 150 一 200 mx572 22。 影像 资料 表明 (网 3a)， 
1994—2008 年 巴 托 拉 冰 舌 继续 后 退 约 410 m, 2008 
年 冰 舌 再 次 前 进 ,至 2021 年 7 月 累计 前 进 约 230 
me。 综 上 所 述 , 近 150 a KEL FEF UK TE BEB 50~60 a 
就 会 发 生 一 轮 前 进 , 但 每 轮 前 进 一 般 只 可 维持 10~ 
20 a ,总体 而 言 , 近 150 a 来 巴 托 拉 冰 川 以 退缩 为 主 
F”, 
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2 巴 托 拉 冰 川 出 水 洞 与 融 水 径流 的 识别 


Fig.2 Identification of drainage system of the Batura Glacier 
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图 3 1972 一 2021 年 巴 托 拉 \ 帕 苏 、 固 尔 金 冰川 进退 变化 
Fig. 3 Terminus positions at different periods for the Batura, Pasu, and Ghulkin Glacier from 1972 to 2021 


帕 苏 冰川 虽 曾 于 1907—1913 年 发 生 过 大 规模 
跃 动 ,导致 冰 舌 前 进 约 1200 m” ,但 此 后 却 长 期 处 
于 退缩 状态 , 仅 在 2000 年 前 后 发 生 过 一 次 前 进 。 
1972 一 1998 年 , 帕 苏 冰 舌 累计 后 退 约 180 m,1998 一 
2002 年 前 进 220 m,2002 一 2021 年 再 次 后 退 约 570 
m。 因 此 近 百 年 来 , 帕 苏 冰川 总 体 处 于 加 速 退缩 状 
态 ( 图 3b)。 

近 百 年 来 , 固 尔 金 冰 川 进退 频繁 , 堪 称 公路 治 
线 最 活跃 的 冰川 之 一 。1885 一 1980 年 固 尔 金 冰川 发 
生 了 3 轮 进退 ,期 间 还 伴随 着 多 次 小 幅度 波动 “。 
1966—1978 ^F. , [5] /K 4 Fa VK TG Ja iB 260 m, 1978— 
1980 EAI X. Bi JE 260 m ,此 后 前 进 速率 减缓 ,1980 一 
1994 年 累计 前 进 100 m'”。 之 后 再 度 豪 退 , 至 2021 
年 累计 退缩 约 390 m。 北 冰 舌 更 加 活跃 (图 3c) , 
1972 一 1990 年 后 退 200 m, 1990 一 1998 年 又 前 进 
170 m。 之 后 , 冰 舌 开始 在 小 范围 内 频繁 波动 ,最 终 
在 1998 一 2018 年 累计 前 进 120 mo 2018 FRYKTE 
再 次 衰退 ,至 2021 年 已 后 退 约 36 mo 


3.2 融 水 通道 变化 

1973 年 6 一 7 月 , 巴 托 拉 河 改道 后 冰川 主 出 水 洞 
由 中 流 线 附近 转移 至 冰 舌 南 侧 , 但 此 后 主 出 水 洞 位 
置 开 始 缓慢 迁移 。2000 年 之 前 迁移 速率 较 慢 ,基本 
稳定 在 冰 舌 南 侧 。2000 年 之 后 ,迁移 速率 明显 加 
快 ,逐步 越过 中 流 线 到 达 冰 舌 北 侧 。 至 2021 年 6 
月 , 主 出 水 洞 已 处 于 冰 舌 西北 侧 , 排 水 道 随 之 延伸 
完全 横 百 于 冰 舌 前 (图 4a)。2021 年 6 一 7 月 , 冰 舌 的 
突然 前 进 挤 压 切 断 了 主 排 水 道 ,引发 了 近 50 a 来 巴 
托 拉 河 的 又 一 次 大 改道 (图 4b)。 由 于 冰 坝 阻挡 , 融 
水 径流 被 迫 在 冰 磺 丘 陵 中 侵蚀 切割 出 一 条 新 的 通 
道 ,直接 导 致 巴 托 拉 河 由 2.9 km 缩短 为 2.5 km, 纵 比 
降 由 65%o 增 大 至 76%o ,因此 流速 显著 增 快 。 

受 地 形 影 响 , 帕 苏 冰川 的 排水 通道 十 分 稳定 。 
但 在 1981 年 之 后 , 冰 舌 附近 在 融 水 补给 下 形成 了 一 
个 小 型 终 磺 阻塞 湖 ,2011 年 实测 面积 约 0.12 km, 
在 水 量 增加 、 冰 舌 崩 塌 、 谷 坡 落石 等 的 共同 作用 下 ， 
阻塞 湖 兽 于 2007 年 7 月 .2008 年 1 月 4 月 等 数 次 发 
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图 4 巴 托 拉 、 固 尔 金 冰川 融 水 径流 排水 通道 变化 


Fig.4 Variations of the meltwater runoff drainage channel of the Batura and Ghulkin Glacier 
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完整 性 已 遭 到 破坏 , 冰 湖 发 育 规模 受到 很 大 限制 。 

固 尔 金 段 是 公路 沿线 受 冰 川 融 水 影响 最 大 的 
路 段 之 一 。 由 于 冰 舌 分 为 南北 两 支 是 融 水 径流 流 
程 极 短 ,不足 1 km, 因 此 冰川 融 水 呈 辐 散 状 , 极 不 稳 
定 且 频繁 发 生 改 道 。 固 尔 金 冰川 主要 排水 通道 约 
有 6 条 ,小 型 通道 则 难以 计数 “ V Henn AL Ke gl 
水 集中 在 M1 通道 , 南 冰 舌 的 M2 通道 也 已 基本 废弃 
(图 4c)。 长 期 以 来 ,M5 为 南 冰 舌 的 主 排水 道 ,但 近 
年 来 M4 的 水 量 也 在 增 大 。M3 曾 在 2014 年 前 后 短 
暂 成 为 主 排水 道 ,但 之 后 水 量 又 迅速 减少 ,M6 水 量 
也 较为 有 限 。 因 此 ,近年 来 南 冰 舌 融 水 通道 以 M5 
与 M4 为 主 ,M3 与 M6 为 辅 。 

3.3 冰 面 高 程 与 流速 

近年 来 “ 喀 喇 昆仑 异常 "现象 备 受 关注 ,但 就 处 

于 喀 喇 昆仑 山西 北边 缘 的 巴 托 拉 地 区 而 言 ,2000 年 
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以 来 的 冰川 物质 亏损 却 有 所 加 速 ”5 。 基 于 KH-9 
DEM 和 SRTM DEM 的 研究 发 现 ,1973 一 1999 年 巴 托 
拉 和 帕 苏 冰川 物质 平衡 分 别 为 0.00+0.10 m w.e. ar 
和 0.05+0.11 m w.e. a' ,表明 2000 年 以 前 冰川 物质 收 
支 基本 保持 稳定 ”'。 但 此 后 的 2000 一 2016 年 , 巴 托 
dir AREA . 固 尔 金 冰 川 表 面 高 程 均 发 生 了 下 降 (图 $a)， 
下 降 速率 分 别 为 0.09+0.29 ma'.、0.15+0.17 m-a Fill 
0.08+0.26 ma ', 尤 以 冰川 末端 表现 最 为 明显 ，。 尽 
管 下 降 速 率 较为 微弱 ,但 这 依旧 表明 2000 年 以 来 巴 
托 拉 地 区 冰川 发 生 了 物质 亏损 。 

长 期 以 来 ,3 条 研究 冰川 均 维持 着 较 高 的 表面 
流速 。1985 一 2018 年 , 巴 托 拉 、 帕 苏 、 固 尔 金 冰川 的 
平均 流速 分 别 为 59.60+0.28 m-a” 106.44+0.44 m-a” 
和 19.50+0.35 ma'。 就 流速 的 空间 分 布 而 言 ( 图 
Sb) ,高 值 区 域 均 出 现在 冰川 上 游 的 粒 雪线 附近 。 
巴 托 拉 冰 川 的 最 大 流速 出 现在 粒 雪 贫 出口 ,第 一 、 
第 二 冰 流 构成 的 冰 瀑 布 处 ,接近 1200 ma!。 帕 苏 
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图 5 巴 托 拉 、 帕 苏 、 固 尔 金 冰川 表面 高 程 变化 (a) 与 年 平均 流速 (b) 


Fig.5 Surface elevation change (a) and annual average flow velocity (b) of the Batura, Pasu, and Ghulkin Glacier 
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冰川 的 最 大 流速 出 现在 相同 的 位 置 , 约 940 ma'。 
固 尔 金 冰川 流速 相对 平缓 ,最 大 流速 也 出 现在 了 冰 
瀑布 以 下 的 附近 位 置 。 相 较 于 巴 托 打 和 固 尔 金 冰 
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跃 动 的 可 能 性 较 低 。 
4.2 融 水 径流 的 影响 因素 
较 高 的 物质 平衡 水 平 导 致 3 条 冰川 的 融 水 径流 


川 , 帕 苏 冰川 流速 较 快 ,意味 着 其 物质 平衡 水 平 更 
高 ,这 是 帕 苏 冰 舌 消融 更 强烈 的 主要 原因 。 


m 


4 讨论 


41 冰川 变化 的 影响 因素 

山地 冰川 是 气候 与 地 形 共同 作用 的 产物 ,但 在 
十 年 或 百年 尺度 上 ,冰川 变化 主要 受气 候 条 件 的 影 
响 25。 冰 川 物质 平衡 直接 由 当年 气候 状态 决定 ,未 
端 进退 却 更 多 地 受气 候 长 周期 波动 的 影响 ,因为 从 
气候 变化 到 冰 舌 波动 需要 一 个 能 量 积聚 和 动力 反 
映 的 过 程 , 即 存在 滞后 性 。 滞 后 期 的 长 短 主要 与 冰 
川 规模 有 关 , 小 规模 冰川 为 3~30 a, 大 规模 冰川 则 会 
ER, 

就 巴 托 拉 、 帕 苏 、 固 尔 金 这 种 大 型 复式 山地 冰 
川 而 言 ,气温 的 长 周期 波动 是 影响 冰 舌 进退 的 关键 
因素 中 。 如 图 5 所 示 ,2000 年 以 来 ,冰川 积累 区 基本 
保持 稳定 ,消融 区 高 程 则 有 所 降低 ,表明 近 几 十 年 
来 罕 萨 河谷 的 升温 直接 导致 了 冰川 物质 亏损 的 加 
FS 。 而 巴 托 拉 、 凰 苏 、 固 尔 金 冰 舌 在 减 薄 的 同时 
发 生 进 退 波动 , 则 充分 体现 了 冰 舌 进退 对 气候 响应 
的 滞后 性 以 及 冰川 规模 对 于 滞后 性 的 影响 。 有 研 
究 表 明 , 近 150 a 来 巴 托 拉 冰 舌 进 退 相 较 于 气温 波动 
存在 60 a 左 右 的 滞后 期 ,冰川 的 数 次 前 进 均 与 气温 
的 阶段 性 变化 存在 很 好 的 对 应 关系 党” 。 基 于 这 一 
判断 ,不 仅 成 功 预测 了 巴 托 拉 冰 舌 在 1970s 一 2000 
年 的 进退 ,还 成 功 预测 了 冰 舌 在 2010s 的 最 新 一 轮 
的 前 进 。 

跃 动 冰川 指 周期 性 发 生 快 速 运 动 的 冰川 , 当 跃 
动 发 生 时 其 速度 往往 是 平时 的 数 倍 至 上 百倍 "。 有 
研究 将 巴 托 拉 第 一 、 第 二 冰 流 等 部 位 的 高 流速 .大 
量 裂 际 与 弧 拱 发 育 视 为 跃 动 的 表现 只 ,但 实际 上 这 


量 明显 高 于 其 他 地 区 的 同等 规模 冰川 , 巴 托 拉 冰 川 
甚至 高 于 近 710 km h REZEKI o EHER 
多 年 平均 流量 为 40.1 m.s…, 瞬 时 最 大 流量 可 超过 
600 ms ,其 中 近 80% 的 径流 量 集 中 在 S—9 H , 近 
60% 集 中 在 7 一 8 月 “2。 巴 托 拉 河 90% 的 水 量 来 源 
于 冰川 融 水 ,因此 径流 量 主 要 受气 温 控 制 , 随 气温 
上 升 而 增加 ,在 夏季 降水 条 件 下 有 时 其 至 会 因 气 温 
下 降 而 减少 ”。 此 外 ,即便 在 7 一 8 月 巴 托 拉 河 的 径 
流量 波动 也 十 分 剧烈 , 常 随 极端 高 温 天 气 的 出 现 而 
陡 涨 引 。 帕 苏 、 固 尔 金 冰 川 融 水 径流 的 年 内 变化 与 
巴 托 拉 河 类 似 , 但 变化 过 程 远 不 及 前 者 剧烈 。 

巴 托 拉 河 径流 年 际 变 化 较为 平缓 ,主要 由 于 : 
(1) 干旱 年 份 降水 偏 少 而 融 水 偏 多 ,湿润 年 份 降水 
偏 多 而 融 水 偏 少 ,起 到 了 天 然 的 调节 作用 ”i,(2) 冰 
舌 区 年 平均 气温 接近 10 % ,导致 消融 十 分 剧烈 , 存 
在 热 融 喀斯特 现象 ,因此 冰 内 裂隙 UC, KREIS) 
清水 能 力也 能 起 到 一 定 的 径流 调节 作用 “2 。(3) 
1970s 以 来 巴 托 拉 流 域 冰川 基本 保持 稳定 ,即便 
2000 年 以 来 物质 亏损 有 所 增加 ,但 对 融 水 径流 的 影 
响 依 旧 有 限 。 实 际 上 不 仅 是 巴 托 拉 河 ,1980s 以 来 
整个 罕 萨 河流 域 的 年 径流 量 也 无 明显 变化 *。 

巴 托 拉 河 在 1973 年 和 2021 年 的 2 次 改道 有 一 
定 的 相似 性 ,改道 均 发 生 在 6 一 7 月 融 水 径流 快速 上 
升 的 时 期 改道 前 主 排水 洞 均 已 发 生长 期 缓慢 的 迁 
移 改道 的 发 生 均 与 冰 舌 突然 前 进 有 关 , 因 此 2 次 改 
道 的 原因 也 是 基本 相同 的 。 如 图 1 所 示 , 巴 托 拉 冰 
川 的 物质 来 源 极 不 对 称 , 冰 舌 北 侧 主 要 由 尤 克 苏区 
效 、 北 \ 西 冰 流 等 构成 . 南 侧 主要 由 第 一 第 二 冰 流 
构成 。 由 于 南 侧 冰 流 拥有 更 大 的 高 差 与 更 开阔 的 
粒 雪 盆 ,因此 物质 供给 量 更 大 ,直接 导致 冰 舌 南 侧 
厚度 与 流速 均 高 于 北 侧 ”“。 冰 川中 下 游 , 南 侧 弧 拱 
与 裸 冰 的 延伸 范围 明显 大 于 北 侧 就 是 二 者 存在 流 


对 于 大 型 山地 冰川 而 言 并 不 罕见 。 对 巴 托 拉 这 种 
巨型 冰川 而 言 , 支 冰 流 跃 动 对 末端 进退 的 影响 也 十 
分 有 限 , 如 克拉 牙 依 拉克 、 费 德 雁 科 等 冰川 支 冰 流 
的 跃 动 并 未 改变 整体 的 退缩 趋势 站 。 帕 苏 冰 川 曾 
于 1907 一 1913 年 , 固 尔 金 冰 川 曾 于 1913 一 1925 年、 
1966—1978 FRÆ EK aJ] ^" , (Hlc 50 a 来 2 条 冰川 的 
积累 区 并 未 发 生 显 著 增 厚 , 因 此 在 近期 发 生 大 规模 


速 差 的 直观 表现 。 在 这 种 流速 差 的 长 期 作用 下 , 主 
排水 洞 受 到 推 挤 而 缓慢 北 移 ,排水 道 随 之 延长 并 横 
吾 于 冰 舌 前 ,最 终 被 突然 前 进 的 冰 舌 切断 而 被 迫 改 
道 。 因 此 , 冰 舌 南北 两 侧 物质 供给 量 的 差异 是 巴 托 拉 
河 改 道 的 根本 原因 , 冰 舌 的 突然 前 进 则 是 直接 原因 ”。 
4.3 对 公路 的 影响 

本 人 研究 在 公路 沿线 选取 了 4 个 稳定 的 参考 点 ， 
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分 别 为 巴 托 拉 桥 、 帕 苏 桥 ,M1 桥 .M4 桥 ,以 参考 点 与 
冰 舌 的 最 近 距 离 来 衡量 公路 与 冰川 的 距离 。 如 图 6 
所 示 , 近 50 a 来 巴 托 拉 冰 舌 与 公路 的 距离 总 体 上 在 
增 大 ,当前 尚 有 1 km 以 上 ,因此 在 短期 内 不 会 对 公 
路 产生 直接 威胁 。 帕 苏 冰 舌 自 2000 年 以 来 不 断 加 
速 后 退 ,与 公路 的 距离 也 在 增 大 ,因此 短期 内 也 不 
A A BUNA S SE TER AE PI AE OE HED 
围 内 波动 ,但 由 于 距离 过 近 , 且 在 将 来 不 能 排除 发 
生 跃 动 的 可 能 ,因此 对 公路 的 直接 威胁 将 长 期 存在 。 
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图 6 1972 一 2021 年 巴 托 拉 、 帕 苏 、 固 尔 金 冰 舌 与 公路 的 距 
离 变 化 
Fig.6 Variations of the distances from glacier tongues of the 
Batura, Pasu, and Ghulkin Glacier to the Karakoram highway 
from 1972 to 2021 


巴 托 拉 冰 川 对 公路 更 现实 的 威胁 在 于 巴 托 拉 
河 新 的 大 改道 。2021 年 6 一 7 月 改道 后 , 巴 托 拉 河 流 
程 缩短 , 比 降 增 大 ,下 切 侵蚀 能 力 增强 , 近 1 a 年 来 ， 
新 河道 明显 拓宽 加 深 (图 7)。 更 关键 的 是 新 河道 下 
游 折 向 南 流 ,与 公路 基本 平行 ,因此 水 流 将 在 转达 
处 直接 冲击 路 基 。 甚 至 完全 可 以 推测 , 正 是 由 于 路 
基 坚 固 才 迫使 水 流 南 流 返回 故道 。 巴 托 拉 河 具有 
流量 大 径流 高 度 集 中 、 陡 涨 陡 落 等 特点 字 ,因此 在 
未 来 的 强烈 消融 季 中 , 随 着 水 量 陡 涨 ,超过 600 m°- 
s' 的 水 流 庄 挟 冰 块 砂 石 奔 腾 下 泄 ,将 对 路 基 产 生 强 
烈 侵 刨 ,对 公路 构成 重大 威胁 。 

帕 苏 冰川 近年 来 消融 强烈 导致 融 水 径流 增 大 ， 
但 由 于 冰 湖 发 育 规模 受 限 ,大 规模 冰 湖 演 决 洪水 对 
于 下 游 公路 的 威胁 趋 于 减弱 。 固 尔 金 冰 川 融 水 轻 
则 在 路 面积 水 结 冰 ,影响 通行 , 重 则 侵蚀 路 基 ,淤积 
桥 底 ,甚至 引发 泥石流 洪水 。 但 在 2008 年 采取 工程 
措施 后 , 固 尔 金 段 的 融 水 灾害 已 得 到 了 有 效 治理 。 
然而 , 固 尔 金 冰川 融 水 径流 的 迁移 摆动 问题 仍 不 容 


74*5l'E 
— 


36°30'45"N 
36?30'45"N 


36*30'N 


74?5l'E 74?52'E 


7 202245 H 27 日 巴 托 拉 冰 舌 区 概况 
Fig.7 Overview of Batura Glacier tongue on May 27, 2022 


忽视 ,一 方面 它 可 能 会 导致 过 水 工程 废弃 ,但 更 重 
要 的 是 径流 量 可 能 会 超过 单一 通道 的 排水 能 力 , 进 
而 引发 新 的 灾害 。 
5 4 it 

中 巴 喀 喇 昆仑 公路 巴 托 拉 段 受 冰川 灾害 威胁 
最 为 严重 ,是 影响 公路 畅通 运行 的 关键 路 段 。 基 于 
野外 考察 与 遥感 影像 等 资料 ,分 析 了 巴 托 拉 、 帕 苏 、 
固 尔 金 冰川 近 百年 来 的 变化 ,变化 原因 及 对 中 巴 公 
路 的 影响 ,主要 结论 如 下 : 

(1) 近 百年 来 ,特别 是 2000 年 以 来 , 巴 托 拉 \ 由 
苏 、 固 尔 金 冰川 总 体 处 于 微弱 退缩 状态 , 冰 舌 发生 
后 退 减 薄 或 基本 保持 稳定 ,对 公路 的 直接 威胁 趋 于 

(2) 在 冰川 消融 加 剧 与 冰 香 波动 的 共同 影响 
下 , 巴 托 拉 段 冰川 融 水 径流 频繁 迁移 改道 ,冰川 融 


水 洪水 . 冰 湖 省 决 洪水 .冰川 泥石流 等 次 生 灾 害 对 
公路 的 威胁 正 日 渐 凸 显 。 


(3) 巴 托 拉 段 冰川 变化 对 公路 最 直接 的 威胁 在 
于 2021 年 6 一 7 月 巴 托 拉 河 新 的 改道 ,改道 后 水 流 
将 直接 冲击 侵蚀 路 基 ,影像 公路 安全 通行 。 
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Recent variations of the Batura, Pasu, and Ghulkin glaciers and their potential 
impact on the Karakoram highway 


LI Zhijie”, WANG Ninglian””, CHANG Jiawen™ 
(1. Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Xi an 710127, Shaanxi, 
China; 2. Institute of Earth Surface System and Hazards, College of Urban and Environmental Sciences, Northwest 
University, Xi’an 710127, Shaanxi, China; 3. Institute of Tibetan Plateau Research, Chinese Academy of Sciences, 


Beijing 100101, China) 


Abstract: The Karakoram highway (KKH) is a strategic channel connecting China and Pakistan and promoting 
the opening up of Xinjiang to the outside world. It is situated on the valley floors of the Pamir- Hindu Kush- 
Karakoram, where the densest mountain glaciers exist in High Mountain Asia (HMA). The KKH has been 
affected by glacier hazards for a long time, especially the Batura- Ghulkin section in the middle reaches of the 
Hunza River. In this study, based on the research documents, field investigation reports, Landsat MSS\TM\ETM+ 
VOLI images, elevation change data sets of HMA, and ITS LIVE glacier surface velocity data, we reconstructed 
the historical changes of the Batura, Pasu, and Ghulkin glaciers, including the glacier area, surface elevation, 
surface velocity, and meltwater runoff over the last 100 years, as well as the impact of glacier variations on KKH. 
The results show that over the past 100 years, the Batura and Pasu glaciers have generally retreated, with the 
distance increased between the glacier tongue and KKH, while the Ghulkin glacier has remained stable. 
Therefore, for the KKH, the current changes in the Batura and Pasu glaciers pose no direct threat. As the Ghulkin 
glacier is too close to the KKH, particularly the northern glacier tongue, which is only about 170 m, the direct 
threat to the highway will persist for a long time. For the Batura glacier, the more realistic threat is the migration 
and reorganization of the meltwater runoff drainage channel (Batura River) from June to July 2021. After the 
diversion, the water flow will cause strong erosion of the roadbed, which must be attended to by the government. 
The Pasu glacier has been retreating strongly in the past decades, and the risk of large-scale glacial lake outburst 
flood has also decreased recently, so the glacier hazard' s threat to the KKH is weakening. With the construction 
of drainage facilities, the meltwater runoff of the Ghulkin Glacier has been effectively channeled, but the disaster 
risk contained in the migration and swing of the meltwater channel can not be ignored. 


Keywords: Karakoram highway; Batura glaciers; meltwater runoff; glacier hazards 


